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Abstract 

The exact solutions of the complete Dirac equation for fermions 
moving in ideal Aharonov-Bohm rings are used for deriving the exact 
expressions of the relativistic partial currents. It is show that as in 
the non-relativistic case these currents can be related to the derivative 
of the fermion energy with respect to the flux parameter. A specific 
relativistic effect is the saturation of the partial currents for high values 
of the total angular momentum. Based on this property, the total 
relativistic persistent current at T = 0 is evaluated giving its analytical 
expression and showing how this depend on the ring parameters. 
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1 Introduction 


The electronic effects in mesoscopic rings were studied by using the non- 
relativistic quantum mechanics n-UDi based on the Schrodinger equation 
with additional terms describing the spin-orbit interaction HU-HS]. 

However, there are nano-systems, as for example the graphenes, where 
several relativistic effects can be observed in the electronic transport. These 
can be satisfactory explained considering the electrons as massless Dirac 
particles moving on honeycomb lattices ini-ra. Other Dirac materials are 
the topological insulators like HgTe and HgTe/CdTe quantum wells with 
low density and high mobility, in which the quantum spin Hall effect can be 
realized [aiE2lE3]. 

Consequently, many studies mm, m-m concentrate on the relativis¬ 
tic effects considering the electrons near the Fermi surface as being described 
by the (l-|-2)-dimensional Dirac equation corresponding to a restricted three- 
dimensional Clifford algebra. However, in this manner one restricts simul¬ 
taneously not only the orbital degrees of freedom but the spin ones too, 
reducing them to those of the S'0(l,2) symmetry. 

Under such circumstances, we believe that there are situations when it is 
convenient to use the complete (l-|-3)-dimensional Dirac equation restricting 
the orbital motion, according to the concrete geometry of the studied system, 
but without affecting the natural spin degrees of freedom described by the 
SL{2, C) group. Thus the polarization effects could be better pointed out. 

Nevertheless, the complete Dirac equation was only occasionally used for 
investigating some special problems of the fermions in external Aharonov- 
Bohm (AB) held as for example the spin effects in perturbation theory [251 
l29[ [30] , the behaviour of the AB fermions in MIT cylinders EH and even the 
AB dynamics using numerical methods. 

The persistent currents in AB quantum rings were recently studied start¬ 
ing with a version of restricted Dirac equation involving a non-Hermitian 
term, introduced by the orbital restrictions [32]. This distort the results 
presented therein as well as other ones based on this approach, even though 
these may outline new realistic effects as in Ref. [33]. For this reason, we 
would like to continue this investigation here, but considering the correctly 
restricted Dirac equation, involving only Hermitian operators, that can be 
obtained easily starting with a suitable restricted Lagrangian theory. 

We discuss this topics showing first that the solutions of the Dirac equa¬ 
tion in the AB rings are determined as common eigenspinors of a complete 
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systems of commuting operators including the energy, total angular momen¬ 
tum and a specific operator analogous to the well-known Dirac spherical 
operator of the relativistic central problems [M]. These solutions can be 
normalized with respect to the relativistic scalar product obtaining thus the 
system of normalized fundamental solutions that allow us to write down the 
exact expressions of the relativistic partial currents and to derive the persis¬ 
tent ones. 

The relativistic partial currents we obtain here are related to the deriva¬ 
tive of the relativistic energies as in the non-relativistic case but, in contrast, 
there appears a crucial difference: in the non-relativistic theory the partial 
currents are proportional to the angular momentum while in our approach 
the relativistic currents tend to saturation in the limit of high total angu¬ 
lar momenta. For this reason we reconsider the problem of the relativistic 
persistent currents at T = 0 proposing an approximative analytical formula 
that matches the numerical calculations with a satisfactory accuracy. 

The paper is organized as follows. In the second section we present the 
relativistic theory of the fermions in AB rings based on a suitable restriction 
of the complete Dirac equation deducing the form of the normalized spinors. 
The next section is devoted to the properties of the partial and persistent 
currents. Finally we briefly present our conclusions. 

2 Dirac fermions in AB rings 

Let us consider a Dirac fermion of mass M moving on a ideal ring of radius R 
whose axis is oriented along the homogeneous and static external magnetic 
held B given by the electromagnetic potentials Aq = 0 and A = Ax . 

The ideal ring is a one-dimensional manifold (without internal structure) 
embedded in the three-dimensional space according to the equations r = R 
and z = 0, written in cylindrical coordinates {t, x) -A {t, r, (f), z) with the 2 ; 
axis oriented along B. Then, it is natural to assume that any held V’ dehned 
on this manifold depends only on the remaining coordinates {t, 0 ) such that 
dri^ = 0 and = 0. These restrictions give the kinetic term, 

5o = y dt , (1) 

of the Dirac action S = So — /3 f dtdcj) '^ 7 '^'^ in the mentioned external mag- 
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netic field, where ip = '^ 17 ° and 

\ 

7'^ = — (—7^ sin0 + 7^ COS0) (2) 

R 

is depending on (p. The notation 13 = ^eBR^ stands for the usual dimension¬ 
less flux parameter (in natural units). 

From this action we obtain the correctly restricted Dirac equation, Eutp = 
Mip, with the new self-adjoint Dirac operator 

Ed = ij^dt + - /3) + ^ > (3) 


whose supplemental last term guarantees that Ed = Ed- This operator 
commutes with the energy operator El = idt and the third component, J 3 = 
Lz+S^i of the total angular momentum, formed by the orbital part L 3 = —id^ 
the spin one >S '3 = ^ diag((T 3 , (T 3 ). Therefore, we have the opportunity to look 
for particular solutions of the form 


'ipEAt,(p) 


N 


y^gi(/>(A+2) 

2 ) 




(4) 


which satisfy the common eigenvalue problems, EDipE,\(t,(p) = MipE,x(t,(p) 
and 

HlpE,x{t, (p) = ElpE,x{t, (p) , J3pJE,x{t, (p) = >^pJE,x{t, (p) , (5) 

laying out the energy E and the angular quantum number A = ±|,±|,... 
whose values are determined by the condition ipE,x(t, (p + 27 i) = ipE,x(t, (p)- 
In this manner we separated the variables remaining with a system of al¬ 
gebraic equations that in the standard representation of the gamma matrices 
(with diagonal 7 °) reads 


/ E-M 0 0 |(A + /9) \ 

0 E-M -j^{X + /3) 0 

0 -;i(A + /3) -E-M 0 

\ ^(A + /3) 0 0 -E-M ) 

This system has non-trivial solutions only for the discrete values of energy 

E^ = ^[M^R^E{l3 + \f]^ , (7) 


( h\ 
h 

9i 
\92j 


= 0 


( 6 ) 
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whose second terms encapsulate the AB effect. For each value Ex we find 
two particular solutions for which 


fi 

/2 




( 8 ) 


where are the usual Pauli spinors of polarization cr = with respect to 
the z axis, 


^1 = 


> = 


(9) 


Thus, we find that for cr = ^ the spinors ® take the form 






2\/'kExR 


0 

0 




( 10 ) 




while for a = — | we obtain the solutions 


UxitA) 


2x/7t ExR 


( 0 \ 

0 


^-iExt 


( 11 ) 


The normalization constants are fixed in accordance to the relativistic scalar 


product 

p27V 

Jn 

( 12 ) 

such that 


{UtUi)=h.y, {(7f,(7J) = 0. 

(13) 


Hence we obtained a pair of fundamental solutions of the same energy and 
total angular momentum but which are not eigenspinors of the operators 
L 3 or S 3 . Therefore, we may ask how these solutions can be defined as 
different eigenspinors of a new operator. The answer is obvious if we observe 
that the desired operator is iP = 2 'y^S 3 which satisfies KU^ = EU^. The 
conclusion is that the spinors are common eigenspinors of the complete 
set of commuting operators {E^, H, K, J 3 }. 
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The operator K introduced above is the analogous of the spherical Dirac 
operator Kd = 7 °(25 ■ L + 1) that concentrates the angular variables of 
the Dirac equation in external helds with central symmetry [M]. Note that 
the genuine three-dimensional operator Kjy cannot be used here because of 
our dimensional reduction such that we must consider the simplihed version 
K The eigenvalues of this operator give the polarization in the non- 
relativistic limit. For this reason we keep this terminology considering that 
the eigenvalues k = ±1 of the operator K dehne the fermion polarization 
with respect to the direction of the magnetic held B. 


3 Relativistic currents 

Using the above results we can calculate the exact relativistic expressions of 
the partial currents on quantum rings, pointing out the difference between 
the genuine relativistic theory and the non-relativistic one. We show that in 
the relativistic approach the partial current tends to saturation for increasing 
A such that the persistent currents at T = 0 will get new properties. 


3.1 Partial currents 


Let us start with the quantum rings where the states of the fermions of energy 
Ex are described by the normalized linear combinations 

^x = c+U^ + c_U^ , |c+P + |c-P = l, (14) 

for which the expectation value of the polarization operator reads, 

{'ipx.K'tpx) = \c+\'^ -\c-\^ ■ (15) 

The partial currents (of given A) coincide in this case with their densities, 
Ix = R'ijjx'y‘^‘ipx, that can be calculated with the help of the matrix (j2). Then, 
observing that 

Uj(t,0)7‘^f/,^(U0) = O, (16) 


we obtain the partial current of a fermion of energy Ex as 

1 (3 + X 1 dEx 


Ix = k+|V+ + |c_|V,- = 


27!R‘ E, 


2tt dl3 


(17) 


^It is known that the forms of such operators depend on the number of space dimensions 
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since = R (t, 0) = R. Thus we find that the partial currents 

are independent on polarization being related to energies in a similar manner 
as in the non-relativistic theory. 

The exact relativistic expressions of the partial currents we obtained here 
depend only on two dimensionless parameters v = (3 + \ and fi = MR (or 
MRc/h in usual units) that are the arguments of the auxiliary function y 
defined as 

This function has the remarkable asymptotic behaviour 

lim x(/i, I/) = ±1, (19) 

u^±oo 


which shows that the relativistic partial currents tend to saturation for large 
values of A. Moreover, for small values of u we can expand 

= -+ 0{iy^). ( 20 ) 

Note that the non-relativistic limit recovers the well-known behaviours 


E\ — M —y E\ — 


u 


2Rii 


h^h = 


1 dEy^ 


271R jj, 277 du 


( 21 ) 


Hereby we conclude that the principal difference is that the relativistic 


partial currents (18) are saturated while in the non-relativistic case we do 


not meet this effect since the function x(/^, ^) is replaced then by the linear 
function that is just its tangent in i/ = 0 as we deduce from Eq. (I20|). This 
result was previously outlined in Ref. [33] but based on the non-Hermitian 
Dirac equation of Ref. [32] • Obviously, the correct saturation effect is given 


by the expression of the partial currents (18) derived here. 


Note that the non-relativistic approximation can be used with a satisfac¬ 
tory accuracy only in the domain where the function x(/r, u) is approaching 

to the linear function -. Our numerical evaluations show that in the domain 

u 

— l/i<i/<l/i the difference |x(/i, z/) — is satisfactory small remaining 
less than 0.05. In addition, we estimate that for \v\ > 5/i the current is 
approaching to its saturation value since |x(/i,±5/i)| = 0.98058. 
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3.2 Relativistic persistent currents 

The above results allow us to derive the total persistent current at T = 0 in 
a semiconductor ring of parameter /i having a even number of electrons Ne 
hxed by the Fermi-Dirac statistics. For the mesoscopic rings with R = lOOnm 
the parameter /i is of the order 10^ — 10^. For example, in a InSh ring of 
this radius, the effective electron mass is M = m* = 0.0135 rUe |35] such that 
/i = 3495. This seems to be the minimal value of /i obtained so far but it 
is possible to obtain smaller values in further experiments with mesoscopic 
rings with R < lOOnm or even with nano-rings having R ~ lOnm. According 
to our estimation, the relativistic effects may be measurable for /i < 10^ 
which means that the actual experiments are approaching to this threshold 
which could be reached soon. 

In all these cases the flux parameter /3 remains very small (less than 
10“®) such that we can neglect the terms of the order 0{/3‘^) of the Taylor 
expansions of our functions that depend on u = X + (3. The total persistent 
current at T = 0 is given by the sum 


I 


E /a= 5: (/a+/-a) 


A=-Af. A=| 


1 


Af 

+/3)+ x{f^,—^ +P)] ( 22 ) 

^=1 


over all the allowed polarizations, A = ±|, ±|,..., where Xp = |(A"e —1). 
Furthermore, by using the expansion 


x(/i, X + f3) + x(h) ~^ + /^) — 2j(/i, A)/3 + 0{(3 ^), 


where 


Hr A) = 




(/i2 + A2)f 

we arrive at the relativistic persistent currents. 


(23) 

(24) 


Xf 


I 5 Imax -jiF ’’ HlH ^ ^ J (h*; '^) ) 


(25) 


A=J 


that can be calculated numerically on computer for any concrete value of /i. 
The function j(/r. A) is simple reaching its maximal value 0.7698 for /i = 
and A = I and decreasing then monotonously to zero when /i and A are 
increasing to infinity. This behaviour is a direct consequence of the saturation 




C 
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0.4474 - 
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c 
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Figure 1: The function c(/r) versus fi calculated for = 0.5/i (a) and 
Xp = when the co-domain is very narrow, 10-5 (b). 


of the partial currents that compensate each other in the saturation zone 
where I\ + I-x —)■ 0. These simple monotony and smoothness properties of 


the function j(/r, A) lead to nice results concerning the values of the sum (25 
when we compute all the allowed contributions. Our numerical examples 
show that when fi is increasing then the functions c(/r) are monotonously 
decreasing tending to an asymptotic value (as in Fig. 1). Consequently, in 
the asymptotic zone, fi > 100, we can use the following approximation 


c(h) - f j{f^,>^)dX= 


giving the dehnitive formula of the relativistic persistent currents 


(26) 


I = 


k 


^TTF 


Xp Ne 

fi 2fi 


(27) 


that reproduces the numerical results with a satisfactory accuracy (under 
10-5). Note that the non-relativistic persistent current, that in our notation 
reads / = klmax, represents a good approximation of Eq. (27) only for 


small values of k (say k < 0.2) for which we can use the approximation 
k{l + /c^)-5 = k + 0{k^) ~ k. 
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4 Concluding remarks 


We outlined here the relativistic theory of the Dirac fermions in ideal AB 
rings. We found the complete system of the commuting operators determin¬ 
ing the fundamental solutions that contains the new operator K which is the 
analogous of the Dirac spherical operator of the central problems. Thus we 
derived the polarized spinors that can be normalized with respect to the rel¬ 
ativistic scalar product. The corresponding relativistic partial currents have 
an interesting behaviour for large values of polarization, tending to the satu¬ 
ration value (27ri?)“^ which depends only on the ring radius. This property 
allowed us to derive the relativistic persistent current at T = 0 giving the 
closed expression (27). 

In our opinion, the form of the relativistic partial current is suitable 
for estimating the total currents even for T > 0 by using the Fermi-Dirac 
statistics. This is because in the saturation domain the currents have very 
small contributions such that the sums become satisfactory convergent and 
can be easily performed on computers or estimated by appropriate integrals 
as in Eq. (26). 

Finally, we observe that the relativistic theory based on the complete 
Dirac equation is able to offer new interesting results in investigating systems 
of low energy as those of the solid state physics that seemed to be destined 
exclusively to the non-relativistic quantum mechanics. 
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